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Abstract 

An analytical-numerical analysis of the singular self-adjoint spectral problem for 
a system of three linear ordinary second-order differential equations defined on the 
entire real exis is presented. This problem comes to existence in the nonlinear field 
theory. The dependence of the differential equations on the spectral parameter is 
nonlinear, which results in a quadratic operator Hermitian pencil. 

1 Introduction. Exact solution to a system of two 
nonlinear wave equations 

The construction of precise regular solutions in systems of interacting classical fields and 
study of their dynamic stability are of great interest in modern nonlinear field theory p]. 

In this paper, we study the problem of stability of such a solution for a system of two 
interacting scalar fields (this solution was reported in paper |2j), the neutral Higgs field 
and a charged linear field (the model was suggested in [3]). In the (H-l)-dimensional 
Minkowski space, the considered field system is described by the Lagrangian 

c = - m\' + w)V2 - (5.0)V2 - h^m' - ^\<f - v'fn. (i.i) 

Here, is a real scalar field; ^ is a complex scalar field; and h, m and v are real positive 
constants. We use the system of units where c = ^ = 1; c is the speed of light in vacuum 
and h is the Planck constant. In this system of units, only the dimension of mass M is 
nontrivial; the dimension of length and time is 1/M ([1], p. 13). In (11.11) 0, ^, and v are 
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dimensionless quantities, and [m] = [h] 
independent variables 

X = {hv/V2) X, i = {hv 
and introduce the new functions 



M . It is convenient to turn to the dimensionless 
2)t (1.2) 



:i.3) 



In what follows, we use variables (11. 2p and functions (II. 3p and omit the tilde over the 
letters. 

The system of the Lagrange- Euler equations for the Lagrangian (II. ip in terms of (II. 2p 
and (II. 3p takes the form 



df^ dx^ 
4 



+ = 0, 



1)0 = 0, t,xGR, 



;i-4) 



:i.5) 



where k is a positive dimensionless parameter and = h"^ /m^. (Here and in what follows, 
the asterisk denotes the Hermitian conjugation.) This system is invariant with respect to 
global (not depending on x and t) transformations ^ — » .^exp(za;) and (j) —0; for the 
motion integrals, it has the energy integral E, charge Q, and topological charge P defined 
by the formulas 



E 



dt 



dx 



+ 



2 \ dt 



+ 



1 (df 

2 \'dx , 



1)^ 



dx. 



oo / 



^ dt ^ dt 



dx. 



:i.6) 



:i.7) 



p = - 

2 



— dx, 
ox 



Thus, the quantities E and Q do not depend on time for those solutions of system (II. 4p . 
(II. 5p . for which the integrals on the right-hand sides of (ll.6p and (ll.7p converge, and P 
(formula (II. 8p ) does not depend on t for any (t){t,x) that has finite limits as oo 
Vt G R. 

By definition, the invariance of Eqs. (I1.4p . (ll.Sp with respect to the transformation 
^ — i> ^exp(iQ;), where a is an arbitrary real number, implies global U(l)-symmetry of 
these equations P], and charge (ll.7p is called a U(l)-charge. 

First of all, we note that system (II. 4p . (II. 5p is known to have the following particular 
solutions: 

(1) the trivial solution ,^ = 0, = 0, which is called a false vacuum since it has a nonzero 
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energy density; traveling wave solutions over the false vacuum ^ = iIj{x ±t), = 0, where 
ip is an arbitrary twice continuously differentiable function (these solutions also have a 
nonzero energy density); 

(2) solutions 

e = 0, 0± = ±1 (1.9) 

with zero energy, which are reffered to as true vacua; and 

(3) solutions of the domain wall type 

^ = 0, (f)^{x) =tanh{V2x/K), (1.10) 

that have the finite energy 

= 4V2/{3k), (1.11) 
zero charge Q, and the nonzero topological charge 

P«, = l 

(of course, the antiwall ^ = 0, 0^ = —(f)w is also a solution to (11.41) . (11. 5p with the same 
energy (11.111) and topological charge = —1). 

Definition 1. A topological soliton (or a domain wait) for system (II. 4p . (II. 5p is a 
solution existing and bounded in the entire space and satisfying the conditions 

lim (j)(t,x) = ±1(= Tl), hm ^(t,x)=0 V t G R, 

i.e., when (l){t,x) has the form of a transition layer between two different vacua such 
that P 7^ (in contrast to a nontopological soliton, when the solution tends to the same 
vacuum value as x — > d=oo, i.e., when 0(t, x) has the form of a splash over the true vacuum 
such that P = 0). 

In these terms, solution (ll.lOp is a topological soliton with zero charge Q. 

Definition 2. A topological soliton is said to additionally bear a ?7(l)-charge if Q 7^ 0. 
Such a soliton is reffered to as a topological Q-ball (in contrast to a nontopological Q-ball 
with P = 0). 

The conditions for the existence of nontopological Q-balls in the Lee-Friedberg-Sirlin 
model [3] and their stability are discussed, along with p], in [1] Chapter 10, and in 
[2]. In particular, it is shown that stable nontopological spherically symmetric and one- 
dimensional Q-balls for such a model exist for large values of the charges, Q > Qc, where 
Qc is a critical charge depending on the parameters of the Lagrangian. However, the 
explicit form of such Q-balls has not yet been found. 

In [2], for system (II. 4p . (II. Sp . the following precise solution of the topological Q-ball 
type was found: 

(pQitjX) = (f)o{x) = tanh(x), (1-12) 
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Ut,x) = J\-^^^^, 0<^'<2, t,xeR. (1.13) 



For this configuration, from fll.6p - (11 .Sp . we obtain 



4/4 N 

i^o-3(^-l), (1.14) 



Qo = 2(|-l), (1.15) 

Po = 1. (1.16) 

When K ^ ^2 the solution ffLT2D . (11331) tends to flTTOD . 

Remark 1. It follows from the above discussions, that the functions <^o = ~4'o 
and ^0 = 'Coexp(ia) are also solutions to system fll.4p . (11. 5p with the same Eq, Qq, and 
Po = -l. 

Remark 2. Taking into account the invariance of Eqs. (I1.4p . (ll.Sp with respect to the 
Lorentz transformations of the independent variables t and x, we find that solution (ll.lOp 
generates for Eq. (11.50 (for ^ = 0) the traveling wave front (l)w{t,x) = tan\^[(^/2 / k) (x ± 



vt)/y/l — f 2] if an initial speed v: < f < 1, is given; solution (11.120 . (I1.13P generates 
for system (II. 4p . (II. 5p the traveling wave of the form 0Q(t, x) = tanh (x ± vt)/^/l — i 



^ / \ /l 1 exp k(t ± vx)/y/l — v^] . , ^ 

t,Q{t,x) = y ^ . with osculations m the ^-component. 

\l K-' 2 cosh[(x ± 'yt)/vl — f^] 

In conclusion of this section, we briefiy discuss physical interpretation (in addition to 
that in [2]) of solution (I1.12p . (I1.13p . This solution possesses properties similar to those 
of possible Q-balls [1] in the (3+l)-dimensional Minkowski space, which are assumed to 
have something to do with the problem of the baryon asymmetry of the Universe. One 
of the mechanisms that can explain the baryon asymmetry has been suggested in [5]. In 
accordance with this mechanism, at the late infiation stages, a condensate is formed that 
can evolve into Q-balls that bear the same baryon charge but are more advantageous 
from the energy standpoint. Note that these Q-balls may continue to exist until now 
and contribute to the dark matter. Moreover, it has been noted in |i3j that such "relict" 
Q-balls may occur crucial factors when studying stability of neutron stars. 

The precise solution of the Q-ball type in the (H-l)-dimensional Minkowski space, 
which was found in [2] and is studied in this paper, is important not only as an approx- 
imation of (3-|-l)-dimensional Q-balls; it is of great interest in connection with studying 
domain walls and processes on them. 
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2 Dynamic stability of the solution and spectral prob- 
lem statement 



From the standpoint of physical apphcations, the problem of the dynamic stability of 
solution fll.l2p . fll.131) is of great importance. A solution is considered to be absolutely 
stable if the energy functional takes its absolute minimum on this solution. For system 
(11.41) . (11. 5p true vacua (11.91) are solutions of this kind. However, these solutions possess 
zero Q- and P-charges. 

Definition 3. A solution ^{t,x)^ <P{tiX) to system (II. 4p . (ll.Sp is said to be absolutely 
stable in a sector {P, Q} if it has the least energy among all solutions with fixed values of 
charges ([II7D and fOj) . 

Solution (1 1.1 op is absolutely stable in the sector {1,0}, since it is the only stationary 
solution of system (II. 4p . (ll.Sp in this sector. 

The question of whether solution (I1.12p . (I1.13P is absolutely stable in the sector 
{Pq,Qq}, where Qq and Pq are defined by (ll.lSp and (I1.16p . respectively, is not easy 
to answer because of the dependence of component (I1.13P on time. 

A suggestive consideration regarding the stability of solution (I1.12p . (ll.lSp might be as 
follows: the parameter affects only the amplitude in (I1.13P : for = 2, solution (I1.12p . 
(ll.lSp turns to the absolutely stable solution (ll.lOp : and, for < < 2, there appears a 
nonzero f/(l)-charge, which usually only stabilizes the solution (see [IJ, Chapter 10). 

The main part of this paper is devoted to studying the dynamic stability of solution 
(I1.12p . (I1.13P with respect to small perturbations (Lyapunov's stability in the framework 
of linear theory). However, first, we present some physical considerations regarding the 
stability of configuration (I1.12p . (ll.lSp from the standpoint of its possible disintegration 
into charged nonlocalized formations. 

1. The problem consists in searching for possible solutions in the sector {Po,Qo} = 
{l,(5o} that are close to solution (11.100 . which is absolutely stable in the sector {1,0}. 
The solutions are sought in the form of small nonlocalized perturbations of the field ^ 
determining the charge Qq. Such solutions, if they exist, may occur equivalent to or more 
advantageous than solution (I1.12p . (I1.13P from the energy considerations. An approximate 



solution to system (II. 4p . (11.50 close to (ll.lOp is sought in the form {5^{t, x), 0^(a;)}, where 
< 1 and 



i.e., only the ^-component in (ll.lOp is perturbed. From (II. 4p . we obtain the following 
equation in 5^: 




(2.1) 



+ 2 



1 - 



1 



0. 



(2.2) 




(2.3) 
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where is a separation parameter and f{x) satisfies the ODE 



ijj 



cosh^ ( \f2 x/ K 



/ = 0, a; G R. 



(2.4) 



For reasons that will be clear from the following discussions, we assume that the solution 
is defined in an arbitrarily large but still finite interval {—L, L). From (11.61) and (11.71) . we 
obtain the following expressions for the charge and energy of solution (12.31) : 



1^ 

Q = 2uj J \f{x)\'^dx, 



(2.5) 



6E 



dx 



+ 2tanh^ (v^x/k) |/p 



dx. 



(2.6) 



In the interval < cj^ < 2, Eq. ([M]) may have only a descrete spectrum. If cu^ is an 
eigenvalue from the discrete spectrum, the corresponding eigenfunction /(x, tu^) belongs 
to L2(— oo, oo) and has the following asymptotics for large \x\: 



-V2 



tu^ \x 



(2.7) 



where A is a normalization constant. Since the perturbation must be normalized by the 
given charge Qo, the localized eigenfunction with asymptotics (12.71) does not generally 
meet the requirement of the perturbation smallness for large Qo- 

The values u"^ > 2 belong to the continuous spectrum. If f{x,uj'^) is an eigenfunction 
corresponding to an eigenvalue belonging to the continuous spectrum, it has the following 
asymptotics for large |x|: 



/(x, uj'^) ~ A cos (^y/uj^ 



\x 



{2.i 



where the phase S is uniquely determined by the boundary condition at x = 0, either 
/(O) = or /'(O) = 0, and the amplitude A is the normalization constant. If normalization 
(12.51) is used, the amplitude A may be selected as small as desired through an appropriate 
choice of L. In addition, /'(x, cu^) ~ for u"^ ~ 2, so that the least value for (12.61) is 
obtained when o;^ ^ 2 in (12. 8p . Since, in the framework of the approximation considered, 
integrals (II. 6p and (ll.7p with the infinite limits diverge on solutions (12. ip . (12.30 . and (12. 8p 
for u"^ > 2, the function /(x, o;^) is set equal to zero outside some arbitrarily large (but 
finite) interval —L<x<L. Note that this expedient is often used when solving physical 
problems (see, for example, [7j, pp. 170-174; pQ, Chapter 10.) 
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Let us multiply fl2.4p by f*{x) and integrate the equation obtained in the interval 

—L < X < L. In so doing, the term -r^f* is integrated by parts with regard to the 

dx^ 

conditions f{^L) =0. As a result, we obtain 

L 



From this and (12.61) . we have 



df 



dx 



+ 2tanh^ (v^a;/^) |/p 



dx = 0. 



6E = 2uj^ J \f\^dx. 



Hence, with regard to (12.51) . it follows that 

6E = uQ > V2Q 



(2.9) 



for cj^ > 2. Finally, we find that the energy of the configuration composed of the domain 
wall (12. ip with charged perturbations (12. 3p over it is the sum of Ey^ and 6E given by 
(II. lip and (12. 9p . respectively; i.e.. 



E = ujQ + 4:V2/{3k)>V2Q + 4:V2/{3k), cu > V2. 



(2.10) 



Note that the energy Eq and charge Qo of the precise solution (I1.12p . (I1.13P are uniquely 
specified by the parameters of Lagrangian (II. ip . At the same time, in the linear approxi- 
mation, configuration (12.10 . (12.30 may have arbitrary E and Q satisfying (I2.1O0 . 

Inequality (I2.10p can be derived in a different way. Since the domain wall (12. ip is 
localized on the interval of length ~ k/v^, to compute the energy of nonlocalized pertur- 
bations of the field ^, we substitute 0^ = 1 and ^ = 6^{t,x) into Eq. (II. 4p . which implies 
the perturbation of vacua (11.90 along the component ^ (here, we rely on the physically 
reasonable assumption that the field ^ is concentrated basically outside the domain where 
the wall is localized). For the function 6^, we obtain the equation (cf. with (12. 2p ) 



9^ _ 



+ 25^ = 0. 



(2.11) 



In order to avoid diverging integrals in the computation of the energy and charge, we 
again consider the perturbations 6^ on a large, but finite, interval —L < x < L. Among 
solutions of Eq. (12. lip corresponding to a given charge Q, the normalized solution 



SC=JQ/{4V2L) exp {tV2t), 



(2.12) 
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has the least energy (note that this solution does not depend on x). Here, the normaliza- 
tion multiplier is determined from the condition that the U(l)-charge for the field 6C, is 
equal to Q, i.e., 

The energy of solution (12.121) is equal to \/2Q. Adding together this quantity and the 
wall energy (II. lip , we obtain the expression coinciding with the right-hand side of (12.101) . 

For the precise solution (11.121) . (11.131) . we obtain from (11.141) and (11.151) the following 
relationship between the charge and energy: 

= 4(Qo + l)/3. (2.13) 

It is easy to see from (12.101) and (12.131) that, for any admissible values of the parameter k, 
the energy of the precise solution is less than the minimal possible energy of the nonlo- 
calized configuration (12.11) . (12. 3p with the same charge. This brings us to the conclusion 
that solution (I1.12p . (11.130 is stable in terms of the disintegration into a nonlocalized 
configuration of the type "wall + plane waves". 

2. The basic objective of this work is to study the dynamic stability of solution (11.121) . 
(I1.13P in the framework of linear perturbation theory. We use the approach similar to that 
employed in studying the stability of localized solutions of some nonlinear wave equations 
in field theory [H] - In particular, it was used in [TU] for studying the stabihty of a 
precise nontopological soliton carrying a U(l)-charge for a complex wave equation with 
fifth-degree nonlinearity. 

Let us set ip = + 5(j) and ^ = + f^^, where 50 and 5^ are small deviations from 
(I1.12p . (I1.13P : note that 6(f)(t,x) is a real- valued function and S^{t,x) is a complex- valued 
function. System (II. 4p . (11.50 reduces to the following linearized system of the differential 
equations in the deviations: 



df^ dx^ 



20^5^ + 400^0^0 = 0, (2.14) 



Qt2 Q^2 + 4eoa'^0 + 40o(eo5r + + ^(30^ - 1)50 = 0. (2.15) 
Similarly [S] - [ID], the solution to Eqs. (12.140 . (I2.15P is sought in the form 

5i = [r]{x) exp(-zAt) + x*{x) exp(zA*t)] exp(zt), (2.16) 

50 = V{x) exp(-iAt) + V*{x) exp{tX*t). (2.17) 

Such a representation of the solution makes it possible to separate variables in (12.140 . 
(I2.15P and to obtain ODEs for the amplitudes of the perturbations rj, x, and V not de- 
pending on t. Indeed, substituting ffLT2ll . (^A^, fl236D . and (12X71) into fICTD and fl27[5D . 
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we obtain the following system of ODEs for rj, x, and V depending on the parameter A: 
,, / ,,o 2\ 11 tanh x 

( cosh^ x) ^ ~^ V^K^ 2 cosh x ' ^ ^ 

\ K"^ cosh X / V ^ 2 cosh x 

Solutions to system (12.181) - fl2.20p are sought in the class of square integrable functions 
defined on the entire real axis satisfying the conditions 

lim rjix) = lim xi^) = li^^ V(x) = 
= lim ri'(x) = lim x'(x) = lim V'(x) = 0. (2-21) 

x^itoo x^±oo x^itoo 

It is required to find the values of A (eigenvalues) for which the singular boundary value 
problem (12.181) - (12.211) has nontrivial solutions (eigenfunctions). By virtue of (12.161) . 
(I2.17p . for any complex eigenvalue A with a nonzero imaginary part, the perturbations 
grow exponentially in time. Hence, the dynamic stability of solution (11.121) . (I1.13P with 
respect to small perturbations of form (12.160 . (12.170 requires that the discrete spectrum 
of problem flZTS]) - fl^D be real. 

In addition, if the ODE system (I2.18P - (I2.20p has a continuous spectrum, it also must 
lie on the real axis of the complex plane A to make solution (11.120 . (I1.13P stable. 



3 Analytic properties of the spectral problem 

Taking into account that the ODE system (I2.18P - (I2.20p is not changed upon the re- 
placement of (x, 7], X, V) by (—X, 1], X, —V) or (— x, —rj, — x, V), we can consider the 
singular boundary value problem on the semiaxis and write it in the following final form: 

,1 ( ,,9 2\ 111 tanh x 

,, / , , 9 2 \ / 1 1 tanh x ^ , , , 

< X < oo, < 2, 
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r/(0) = x(0) = V^'(O) = or ry'(O) = x'(0) = V'(O) = 0, (3.4) 
lim riix) = lim n'ix) = 0, (3.5) 

lim y(x) = lim y'(x) = 0, (3.6) 

lim F(x) = lim V'(x) = 0. (3.7) 

x^oo z— >oo 

This problem contains terms depending nonlinearly on the spectral parameter A and, 
generally, may have complex eigenvalues. We seek for nontrivial solutions to this problem 
in the class of complex-valued functions satisfying the condition 

[r]*{x)ri{x) + x*{x)x{x) + V*{x)V{x)] dx < oo. (3.8) 

Note that the singular boundary value problem (13. ip - (13. 7p is correctly defined in 
terms of the number of the boundary conditions for large x if and only if the singular 
Cauchy problem (13. ip - (13. 3p . (13. 5p - (13. 7p admits a three-parameter family of solutions 
at infinity. Since system (13. ip - (13. 3p is asymptotically equivalent to a system with 
constant coefficients, each of the decoupled second-order ODE for large x must have 
a one-parameter family of solutions vanishing at infinity. This results in the following 
requirements on the location of the eigenvalues A in the complex plane: 

1 + 2A- A^ ^ (-00,0], 1-2A- A^ ^ (-00,0], 8//€2 - A^ ^ (-oo, 0], (3.9) 

i.e., they do not lie on the nonpositive real semiaxis. Then, for sufficiently large x, x 3> 1, 
the three-dimensional linear subspace generated in the phase space of system (13. ip - 
(13. 3p by the values of the solutions of the singular Cauchy problem (13. ip - (13. 3p . (13. 5p - 
(13. 7p is given in the form 

7]'{x) ^ -V1 + 2A- A2r7(x), (3.10) 



X'{x) ^ -Vl-2A-A2x(x), (3.11) 

V'{x) ^ -^/s/^^^^Vix), (3.12) 

where the roots are assumed to have positive real parts (see pQj for detail). 
Hence, the following assertions are valid: 

(1) real eigenvalues of problem (13. ip - (13. 8p . if exist, satisfy (by virtue of (13. 9p ) the 
inequalities 

- V2 + 1 < A < V2- 1; (3.13) 

(2) the geometric multiplicity of any eigenvalue of problem (13. ip - (13. 8p cannot be 
greater than three; i.e., not more than three linearly independent eigenfunctions may 
correspond to one eigenvalue; 
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(3) system fl3.ip - fl3.3l) has a continuous spectrum lying on the real axis of the complex 
plane A in the intervals (—00, —^J2 + 1) and (v^ — 1, 06) (conditions fl3.10l) . fl3.11l) . and 
f l3.12p turn to conditions of the radiation type for A G (—00, — -\/2+l) and A G (l + -\/2, 00), 
for A G [—00, — 1) and A G (-\/2 — 1, 00), and for A^ > S/ > 4, respectively); these 
real intervals of the continuous spectrum do not affect stability of solution (11.121) . (I1.13P 
with respect to small perturbations of form (I2.16p . (12.170 . 

Note also that, if A is a purely imaginary eigenvalue of problem (13. ip - (13. 8p . then 
?7 = X*5 and V is a. real- valued function. 

These assertions can be supplemented by the following facts and estimates. 

1. Case of K^=2. For = 2, system (13.10 - (13. 3p is decoupled into the three 
second-order ODEs: 

7]" + [2/ cosh^ X - (1 + 2A - A^)] r] = 0, (3.14) 



2/ cosh^ X - (1 - 2A - A^) x = 0, 



V"+ 6/ cosh^ a; - (4 - A^ 



V = 0, 



(3.15) 
(3.16) 



each of which can be reduced to a hypergeometric equation. Indeed, consider the equation 



+ 



s{s + l) 
cosh^ X 



— e 



^/^ = 0, —00 < X < 00, 



(3.17) 



and introduce the change of variable y = tanhx. For the function ip{y), we obtain the 
equation 



d_ 

dy 



1 ^\^^ 



+ 



Sis 



ij = 0, 



-Ky <1, 



which reduces to a hypergeometric equation by means of the substitution 4j{y) = {1 — y 
and the change of variable (1 — y)/2 = z. Assuming that w = w{z), we obtain 

z{l - z)w" + [{e + 1) - {{e - s) + {e + s + 1) + l)z]w' - {e - s){e + s + l)w = 0, 

<z <1. 



This equation has the two independent integrals [12] : 

wi{z) = F{e — s, e + s + 1, e + 1, z), W2{z) = z^'^F{s + 1, — s, 1 — z) 
where F is the hypergeometric function 

r(or + A;) 



Via] 
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Hence, returning to fl3.17p . we obtain 



i'lix) = — 77— F (e — s, e + s + 1, e + 1, (1 — tanhx)/2) , 
cosh X 

ip2{x) = 2'^ exp (ex) F {s + 1, —s, I — e, (1 — tanhx)/2) . 

Then, only ipi belongs to L2(— 00, 00) under the condition that e > and e — s = —n, 
where n = 0, 1, 2, . . ., when F becomes a polynomial of the n-th degree in tanhx and, 
hence, has finite limits as x —>■ ±00. 

For example, for Eq. (13.141) . we have s = 1, and e = vT+2A — = 1 — n > 0. 
Hence, it follows that = 0, Ai = 0, A2 = 2, and ipi{x,Xi) = ipi{x,X2) = ri{x) = 
1/ coshx. Equations (13.151) and (I3.16P are treated in a similar way. Hence, in the case of 

= 2, all points of the discrete spectrum of problem (13.11) - (13.81) and all eigenf unctions 
corresponding to them can be found (all eigenfunctions are accurate up to the normalizing 
multipliers): 

Ai=0: r] = 0, x = 0, 1/ = l/cosh2x; (3.18) 
Ai = 0andA2 = 2: r7 = l/coshx, x = 0, V = 0; (3.19) 
Ai = and A3 = -2 : rj = 0, x='^/coshx, V = 0; (3.20) 

A4,5 = ±-\/3 : 1] = 0, X = 0, = tanhx/ coshx. 

Thus, for = 2, the eigenvalue A = Ai = of problem (13.11) - (13.81) has the maximal 
possible geometric multiplicity equal to three, whereas the multiplicity of the eigenvalues 
A25 A3, A4, and A5 is equal to one. 

2. Case of 0\k,^\2. It is easy to check that, for any k satisfying the inequalities 
< < 2, A = Ai = remains an eigenvalue of problem (13.11) - (13.81) . and there are, at 
least, two eigenfunctions corresponding to it; namely, 

Ai = 0: r/ = -X = 1/coshx, V = 0; (3.21) 



/ 1 1 tanh X 1 , ^ 

Note further that, if (77(0;), ^(a;). A) is a solution to problem (13. ip - (13. 8p . then 
{r]*{x),x*{x),V*{x), X*), (x(a;), ?7(a;), K(x), -A) and {x*{x),V*{x)yV*{x), -X*) are also its 
solutions, so that it is sufficient to consider, for example, the following range of A: 

Re A > 0, Im A > 0. (3.23) 

Like in [H] - [ID], we write system (13. ip - (13. 3p in the form 

A^/^ - 2AD^ - i/^ = 0, < X < 00, (3.24) 
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where \l/ = (77, x, V)'^ (T denotes the transposition) and /, H, and D are operators defined 
in the space of three-component complex- valued twice differentiable functions bounded 
on the semiaxis and satisfying conditions (13 ■4p - (13. Sp . \E': [0, 00) — > C^, with the scalar 
product defined as 



^*(a;)*(x) dx. 



Here, / is the identity operator and D and H have the form 





1 










Hn 





Hi3 




D = 





-1 





H = 









(3.25) 















H31 


H32 


H33 





where 



H- 



H22 



33 



dx"^ 



cosh X 
1 



1, 



+ 



cosh^ X k'^ 



(3.26) 



H 



13 



31 



H' 



23 



32 



/ 1 1 tanhx 



2 coshx 

Taking the scalar product of (I3.24p with \E' and assuming that \E' is normalized such that 
(^, ^) = 1), we obtain 

A2-2A(^,L>^)-(^,i/^) =0. (3.27) 

H and -D are Hermitian operators. The operator D is obviously Hermitian; H is Hermitian 
because the operator —d'^/dx^ is Hermitian. The latter follows from the fact that the 



dm 

substitutions W* — — 

dx 



and ^— \1' 

ax 



turn to zero on the functions satisfying conditions 



dMD - ([32D. Then, (^,L>^) and ($,i^$) are real numbers. From it follows that 

A = (^, D^) ± v/(^,£'^)^ + (^,//^). (3.28) 

From this equation, we can obtain estimates for an eigenvalue with a nonzero imaginary 
part. Indeed, in this case, it follows from (13.280 that 



ReA=(^,D^), ImA = V-(^,^^) - (*,^*), |A| 



(^,i/^). (3.29) 



We have the following estimates for the real and imaginary parts of A: 



iReAl = \(m.Dm) 




00 



< 



dx {r]*r] - x*x) 
< j dx \7]*r] - x*x\ < j dx {r]*r] + x*x) < 1; 



(3.30) 
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(3.31) 



where /imm is the least eigenvalue of the operator H. Let us obtain an estimate for fimin- 
From (IsSD, dSSHD, we have 

oo 


oo 

+ Jdx [7]*H,,V + x*H23V + V*HsiV + V*H,2X] ■ 



We estimate for the first integral is obtained by dropping the positive terms: 



dx 



cosh^ x) ^ ^ 
+ V* -— + 



rf2 ^_ 2 



dx"^ k"^ \ k"^ 



+ 2 



cosh^ X / 



V 



cosh X 



dx 



+ 2 tanh^ x - l\ri + x \ -— + 2 tanh^ x - 1 x + 



dx"^ 



+ 1/* 



c/x^ 

— + — + 2)tanh2x-2 1 
dx^ V 



> 



>- J dx [ri*ri + x*X + 2l^V] > -2 J dx [ri*ri + x*X + V*V] = -2 ^) = -2. (3.32) 



tanhx 



Taking into account that 



coshx 



< -, we obtain the estimate for the second integral 



dx4dX-l^-^[V*{x + v) + {x + vrV] 



k"^ 2 coshx 
tanhx 



coshx 



[v*ix + v) + {x + vrv] 



< 



< 



<^\l^-lJdx\V*{x + v) + {x + vrV\. 



(3.33) 



Further, applying the transformation 

r] = {u + w)/2-t/V2, x= {u + w)/2 + t/V2, V = {u-w)/V2 



(3.34) 
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and noting that it is unitary (does not change the value of the scalar product (\E', \E'), we 
obtain 



oo 

\u*u — w*w\ < 



J dx \V*{x + V) + {X + TlTV\ = V2jdx 



oo oo 

<V2jdx {u*u + w*w) <V2jdx {u*u + w*w + t*t) = V2 ^) = V2. (3.35) 





From f l3.33p and (13.351) . we obtain 





From this inequality and fl3.32p . it follows that 



< 2 J2/k^ - 1. 



^^imin >-2(l + - 1 j . (3.36) 

Using (I3.29P - (13.310 and (I3.36p . we finally obtain the following estimates: 



I A| < V^JI~. <\2il + J2/k^ - 1 , (3.37) 



|ReA|<l, |Im A| < V-I^min < y 2 (^1 + - ij . (3.38) 

4 Numerical study of stability 

4.1 Refiniment of the localization domain for the eigenvalues of 
problem - (^M 

Along with analytic estimates (I3.37p . (13.380 . the localization domain for the eigenvalues 
of problem (13.10 - (13. 8p can be determined numerically starting from the estimate 



|A| < \/-/imm(/«)- (4.1) 



To find the numbers /imm(^); we consider an accompanying eigenvalue problem for the 
operator H (see (13.250 . (I3.26P ). i.e., the yU-parameterized ODE system 

jjjxif -H^ = 0, < s < oo, (4.2) 
subject to the boundary conditions (13.40 - (13.70 and requirement (13. 8p . 
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In order to solve this problem numerically, it should be reduced to a problem defined 
on a finite interval. As noted in Section 3, the condition that the solutions belong to the 
space L2[0, oo) is equivalent to the requirement that the values of the solutions for large 
X belong to a three-dimensional subspace of C^, which, up to exponentially small terms, 
is defined as 

v'ixoo) = -\Jl - /i ^?(a;oo), (4.3) 
X'(xoo) = -V^l-/i x(a;oo), (4.4) 
V\x^) = Vix^), (4.5) 

where /i < 1 and > 1 (cf. with flTTUj) - flXT^ l 

Problem (14. 2p . (13. 4p . (I4.3p - (I4.5I) is self-adjoint. Taking into account item 1 of Section 
3, we immediately find that, for k = a/2, /i = is its only eigenvalue (of multiplicity 3), 
and that the corresponding eigenfunctions are given by (I3.18P - (I3.20p . For k different 
from \/2, the eigenvalue /i = has multiplicity 2, with the corresponding eigenfunctions 
being given by (I3.2ip and (I3.22p . and, as computations show, the third zero eigenvalue 
starts to move to the region of negative values. The eigenvalue is sought numerically by 
transferring boundary conditions (14. 3 p - (14. 5 p from the point x = Xoo to the point a; = 
and examining the changes of sign of the resulting determinant. Note that, in this case, 
a variant of the sweep method is used [13] (with regard to the studies on the robust use 
of orthogonal sweep methods in singular problems [Ti]). 

Our computations show that, when k varies from ^/2 to zero, the double zero eigenvalue 
remains unchanged, whereas the single eigenvalue (which is just the eigenvalue fimin we 
are interested in) moves to the region of negative values (Fig. 1). The comparison of these 
results with (I3.36P shows that the bound obtained analytically is too high. For example, 
for K = 1.2, the value of y/—^min found numerically is equal to 0.45, whereas (I3.36P 
yields 1.8. Nevertheless, for k ^ 1.23, the magnitude of ^—[imin is greater than -\/2 — 1 
(see estimate (13.130 and Fig. 2), which complicates the localization of the eigenvalues of 
problem (I3.ip -( !378l) (see below). 



4.2 Search for eigenvalues of problem ( 13.11) — (13.81) in the sector 

03]), (I323D 

Now, let us consider problem (13.10 - (13. 8p itself. Its reduction to a problem defined on a 
finite interval is similar to that for the previous problem for the ODE (14. 2p . Representation 
(I3.10p - (13.120 yields the following boundary conditions at x = x^o'- 



X (-^oo) 



-Vl + 2A-A2r7(xoo), 



Vl-2A-A2x(a:oo), 



(4.6) 
(4.7) 
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Figure 1: the dependence of Umin from k. 



(4J 



where the roots are assumed to have positive real parts. 

The nonsingular spectral problem (13 .ip - (13 ■4p , (14 .6 1) - (14 .Sp was studied numerically by 
a method based on the generalization of the argument principle (methods of localization 
of descrete spectrum points based on the argument principle and its modifications are dis- 
cussed, for example, in [15] - [18] and in references therein) with the use of the differential 
sweep [I3] (as noted earlier, the problems of the robust application of the method from 
|13j . as well as other modifications of the sweep method, to singular eigenvalue problems 
are discussed in [H]). Conditions (14.61) - (14. 8 p are represented in the form 



V?(xoo, A) W{xoc,\) = 0, 



(4.9) 



where W = {rj, i]', x, x'? V, V')'^ and 

(p{Xoo, A) = 











VI + 2A - A2 1 

VI - 2A - A2 1 







- A2 1 



(4.10) 
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Figure 2: the dependence of y/—iJmin from k. 



As a result of the sweep, we obtain a condition at the point a; = with the 3x6 
matrix ip{0,X) equivalent to fl4.9p . fl4.1Up . If we augment thi matrix up to a square one 
by adding the lower matrix block corresponding to conditions (13.41) . eigenvalue problem 
can be formulated as follows: 

a number A is an eigenvalue <^==^ the determinant of the resulting matrix is equal to 
zero; 

the geometric multiplicity of an eigenvalue (i.e., the number of linear independent eigen- 
functions) is equal to the defect of the resulting matrix. 

Note that the coefficients of the matrix ip{0, A) and, hence, the determinant of the 
resulting matrix, are not analytic functions of A. Nevertheless, as shown in [16], the 
argument principle can still be used for determining the number of the eigenvalues in a 
given region of the spectral parameter A. The search can rely on the determinants M245 
for the first condition and M136 for the second condition in (13. 4p . respectively, where 
Mijk is the minor of the matrix (f{0, A) composed of the columns i, j, and k. The zeros 
of det(Mjjfc) correspond to the eigenvalues of problem ( 13.11) - (13. 4p . (14. 6 p - (14. 8p . The 
algebraic multiplicity of an eigenvalue is the multiplicity of the corresponding zero of the 
determinant. 
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From the computational results, it follows that when k is close to v^, the search 
domain for unstable eigenvalues is a circle of radius a/— /imm centered at the origin (see 
(14.11) ). When k approaches zero, the search domain is the same circle with the cuts along 
the real axis from \/2 — 1 to the right and from —\/2 + 1 to the left (see (I3.13P ). 

The critical difficulty associated with this problem is that the points lying on the cuts 
belong to the continuous spectrum (when A is close to the cuts, the computation becomes 
unstable). Thus, the methods of the localization of the eigenvalues described in [T3]- 
[18] cannot directly be applied to our problem, since the resulting determinant (or some 
other function used in the argument principle) cannot continuously be extended to the 
boundary of the domain. (Note only that, if there were descrete eigenvalues appearing 
in the intervals of the continuous spectrum, they would develop somehow when k varies; 
however, the computation did not reveal anything of this kind near the cuts when k 
varied.) 

From the problem symmetry, it follows that the eigenvalues with nonzero real and 
imaginary parts are located in the complex A-plane by quadruplets (A, —A, A*, —A*) and 
those with only real or imaginary nonzero parts, by pairs on the real or imaginary axis, 
respectively. 

To study perturbations of the degenerate eigenvalue A = and reveal the presence / 
absence of complex eigenvalues, we made computations on the following contours in the 
A-plane: 

circles of radius varying from e to \/2 — 1 — £, where e is a small number, centered at 
A = 0; 

circular sectors of the angle 90° with the vertex aX X = is and radius ^/—iJimin located 
in the first quadrant (the eigenvalue symmetry property is used); 
intervals lying on the real or imaginary axis; 
separate contours outside the real axis. 

The computations show that, in the degenerate case of k = a/2, the algebraic mul- 
tiplicity of the eigenvalue A = is equal to four, and its geometric multiplicity is equal 
to three (cf. with (13.181) - (13.201) ). The difference in the two multiplicities is explained 
not by the presence of Jordan block (which is impossible in view of the self-adjoint ness 
of the boundary value problem) but by the fact that Eq. (13. 3p contains the square of the 
spectral parameter A. 

When K varies from \/2 to 0, the following picture is observed. As before, the alge- 
braic multiplicity of the eigenvalue A = is equal to 4. Thus, we may conclude that the 
eigenvalue A = is not split into several eigenvalues. Moreover, when k is close to a/2, 
it follows from the results obtained that the problem has no other eigenvalues. When we 
move from n = a/2, the geometric multiplicity reduces from 3 to 2 (instead of eigenfunc- 
tions (13.191) . (I3.20p we have (I3.2ip ). The reduction of the geometric multiplicity upon the 
perturbation of a parameter is a typical phenomenon for quadratic pencils (in contrast to 
linear ones). Eigenf unctions fl338D - I^M) for k = ^2 and (K2T\\ . (Km for other k are 
found numerically with accuracy up to 10^^. As n approaches zero, problem (13. ip - (13. 8p 
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becomes more stiff. In tlie computations with fixed (double) relative accuracy, the lowest 
value of K was as small as 0.05. 

For additional stability control for the eigenvalue A = 0, we used the following con- 
sideration. From the symmetry of the eigenvalues and the existence of the eigenfunctions 
f l3.2ip . f l3.22p . it follows that, at least, the property of being repeated for the eigenvalue 
A = is conserved. Hence, the perturbation could accur only through splitting of two 
symmetric eigenvalues along one of the axis. We have checked numerically that this effect 
did not take place (with regard to estimates fl3.13p and f l3.38p ). No other eigenvalues have 
been found in the numerical experiments. 

It may be concluded (with the greatest certainty for the interval k: 1.23 ^ k < \/2) 
that the singular boudary value problem, (13. ip - (13. 8p has only one eigenvalue A = of 
algebraic multiplicity 4 and geometric multiplicity 2 with the corresponding eigenfunctions 
(I3.2ip . (I3.22p . No complex eigenvalues have been found in the admissible domain (14.10 . 

5 Conclusions 

The analytical and numerical studies reported in this paper allow us to conclude that 
solution (I1.12p . (I1.13P is dynamically stable with respect to small perturbations of form 
(I2.16p . (I2.17P for < /t < a/2. This statement has the greatest certainty in the interval 
k: 1.23 ^ K < Moreover, solution (I1.12p . (I1.13P in the course of its evolution cannot 
split into small nonlocalized oscillations (along the component ^) over kink (11.1 OP since it 
is not advantageous from the energy standpoint. The question of whether the topological 
Q-ball (I1.12p . (I1.13P is absolutely stable in the sector {P^Qo} is still open. 
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